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We investigate the curvature-dependence of water dynamics in the vicinity of hydrophobic spheri¬ 
cal solutes using molecular dynamics simulations. For both, the lateral and perpendicular diffusivity 
as well as for H-bond kinetics of water in the first hydration shell, we find a non-monotonic solute- 
size dependence, exhibiting extrema close to the well-known structural crossover length scale for 
hydrophobic hydration. Additionally, we find an apparently anomalous diffusion for water moving 
parallel to the surface of small solutes, which, however, can be explained by topology effects. The 
intimate connection between solute curvature, water structure and dynamics has implications for 
our understanding of hydration dynamics at heterogeneous biomolecular surfaces. 


One of the greatest advances in our understanding 
of the hydrophobic effect is the recognition that the 
hydration structure and thermodynamics of apolar so¬ 
lutes is qualitatively length scale dependent [1-11]. The 
microscopic reason is that water structures very differ¬ 
ently at small (convex) solutes, where the bulk H-bond 
network is only moderately deformed, as compared to 
large solutes, which significantly distort the tetrahedral 
bulk structure. The structural crossover happens at sub¬ 
nanometer length scales and has important implications 
for the interpretation of the structure and thermodynam¬ 
ics of hydrophobically-driven assembly processes [1, 12], 
such as protein folding and association [13-15]. 

The dynamics of the hydration layer that surrounds 
molecular self-assemblies and proteins in solution has at¬ 
tracted plentiful interest in the last decade [16]. Solute 
fluctuations and hydration dynamics are understood to 
be highly coupled with important consequences to bio¬ 
logical function, such as enzyme catalysis and molecular 
recognition in binding [16-26]. Despite the obvious im¬ 
portance of the solute chemical composition, apparently 
the intrinsic topological and geometric features play an 
important role as well [19-21, 27], possibly even leading 
to anomalous diffusion behavior [20]. Therefore, and due 
to the established fact that water considerably restruc¬ 
tures at radii of curvature close to the sub-nanometer 
scale, a natural question to ask is, how does the water 
structural crossover affect the dynamics of the hydration 
layers in the solute vicinity? 

One of the first simulation studies of curvature effects 
on hydration dynamics was performed by Chau et al. for 
three solute radii between 0.35 and 0.8 nm [28]. They 
found a slowdown of the diffusion of water in the first 
hydration shell relative to the bulk with an apparent min¬ 
imum for the intermediate solute size. This interesting 
finding was not commented on, probably due to the little 
amount of data and statistical uncertainty of the results. 
Further, a slowdown of water reorientational dynamics 
was found compared to bulk, an effect that decreased 
with solute size [28]. The reorientational slowing-down 


has been recently explained by excluded-volume effects 
on the H-bond exchange dynamics [29, 30]. In disagree¬ 
ment with the simulations with Chau et al, however, the 
excluded-volume concept predicts a monotonic increase 
of the reorientation times [29], an apparent controversy 
which has not been addressed in literature, yet. 

In this letter, we report on systematic molecular dy¬ 
namics (MD) simulations of water around hydropho¬ 
bic spherical model solutes of varying radius between 
0.3 and 2.1 nm and investigate intrinsic curvature ef¬ 
fects on water hydration dynamics. We find that 
the perpendicular and lateral diffusivity of hydration 
water exhibit a non-monotonic curvature dependence, 
with temperature-dependent minima located close to the 
structural crossover length scale [7, 31]. Both, curva¬ 
ture and temperature dependence, are strikingly similar 
to that of entropy during hydrophobic solvation [31] and 
thus allows strong conjecture that the hydrophobic effect 
extends beyond thermodynamic into dynamic anomalies. 
Furthermore, we find that the intriguing curvature de¬ 
pendence of the diffusivity is found to be related to a 
non-monotonic dependence of H-bond life times on cur¬ 
vature, exhibiting a maxima at the crossover length scale. 
This finding reconciles previous, apparently contradict¬ 
ing results from simulations [28] and theory [29, 30] 
on reorientation times, as we find both predicted be¬ 
haviors but in different length scale regimes. The in¬ 
timate connection between solute curvature and water 
structure and dynamics should be of fundamental im¬ 
portance for biomolecular function mediated by hetero¬ 
geneous biomolecular surfaces [16, 17, 19-21, 27]. 

Our MD simulations are performed each containing a 
single and fixed hydrophobic model solute solvated in 
SPC/E water [32] using the Gromacs simulation pack¬ 
age [33]. The solute-water interaction is mediated by a 
shifted Lennard-Jones potential Ulj[t') = 4e[(a/r') 12 — 
(cr/r') 6 ], whereas r' = r — r 0 describes a coordinate shift 
by r 0 . With such a shifted potential, it is guaranteed 
that by changing the solute size only intrinsic curvature 
effects are probed, not, e.g., those of an additionally vary- 
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Fig. 1: Left: MD simulation snapshot of a hydrophobic sphere 
with shift radius ro = 8 A in explicit water. Right: Radial 
distribution function (RDF) g so (r) of water around differently 
sized hydrophobic model solutes. 

ing dispersion attraction. Values of er = 3.768 A and 
e = 1.197 kJ/mol are chosen from a model for methane 
to mimic a reasonable dispersion attraction with a well- 
defined first solvation shell [6]. We choose shift radii r 0 = 

0 A, 1 A, 2 A, 4 A, 6 A, 8 A, 10 A, 12.5 A, 15 A and 17.5 A in 
ten separate simulations. As a limiting case of a hy¬ 
drophobic surface with zero curvature, the box delimit¬ 
ing walls in the x — y-plane in a pseudo-2D simulation are 
chosen to interact with water by a 12-6 potential in z- 
direction Ui 2 -e( z ) = 4e[(cr/^r) 12 — (cr/ 2 :) 6 ]. After lOOps of 
Gibbs ensemble ( NPT ) equilibration, canonical ( NVT ) 
production runs are performed up to 200 ns at a tem¬ 
perature of T = 300 K and with N = 6000 to 12000 
water molecules. Further details on the simulation tech¬ 
nique can be found in the Supporting Information (SI). 
A simulation snapshot and water oxygen density profiles 
around all the solutes are presented in Fig. 1. 

We first characterize the water diffusion parallel to the 
solute surface. For this, we calculate the lateral mean 
square displacement (MSD) of the arc length a water 
oxygen has traveled in time t, via 

(S(t) 2 ) = Rl s -([6(t' + t)-9(t')] 2 ) Rl . (1) 

Here R avg = f^ 1 d rr ■ g S o{f) denotes the average dis¬ 
tance of the water molecules inside the first hydration 
layer, which is delimited by the location R\ of the first 
minimum in the RDF (cf. Fig.l), to the solute center. 
The time average over starting times t' is taken 

only in the first hydration shell. The variable 9 it) is 
the azimuthal angle of the water oxygen-solute distance 
vector at time t. For water molecules diffusing at a 
wall, Eq.(l) reduces to the usual two dimensional MSD 
([Ax(t) + A y{t)] 2 ) r 1 in the x — y-plane. 

The lateral MSDs for water around all solutes are 
shown in Fig. 2(a) along with those of water at the planar 
hydrophobic surface and in bulk. Comparison to the bulk 
MSD indicates a crossover from below to above bulk wa¬ 
ter self-diffusion with increasing slopes with growing so- 



Fig. 2: (a) MSD of water molecules parallel to differently 

curved hydrophobic surfaces and MSD of water in bulk (thick 
black dashed line), (b) The data from the top panel (a) di¬ 
vided by At. The Inset shows the curves for water at the 
five smallest solutes along with the corresponding fits from 
Eq. (3.) 

lute size. Additionally, the hydrating water of the small¬ 
est solutes exhibits non-linear, anomalous behavior in t. 
This can be better recognized by the negative slope for 
small solutes in Fig. 2(b) where we plot the MSD divided 
by time, (S(t))/4t. Linearity of the MSDs is restored 
for larger solutes, converging towards the limiting MSD 
r 0 —> 00 of the planar hydrophobic surface. 

The apparently anomalous diffusion behavior for small 
solutes can be explained by intrinsic curvature effects, 
which modify the standard 2D diffusion law [34]. Here, 
the probability distribution function (PDF) of diffusion 
on spherical surfaces reads 

P{9,t) = ^ J g s i n (fl) . fv (2) 

r 

where r = 2D^t/R 2 vg , and N{t) is a normalization con¬ 
stant. By expanding Eq. (2) up to second order for small 
6 1 we find that the second moment of the PDF can be 
written for small solutes as 

4 Dlt 2 

(6 2 )Rl vg K4D llt ---^. (3) 

On surfaces with high curvature (R 2 vg D||t) the sec¬ 
ond term on the r.h.s. of Eq. (3) slows down the MSD 
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Fig. 3: Parallel and perpendicular diffusivities D\\ (green 
dashed lines) and D± (red solid lines) scaled by bulk water 
diffusion Dbuik against solute size Rg for T — 300 K. For the 
four smallest cavities the temperature effect on D± and D\\ 
is shown for T = 260 K, 280 K and 320 K. The grey dashed 
horizontal line represents bulk water diffusivity. The values 
corresponding to diffusivity at the limiting case of zero cur¬ 
vature (Rg —» oo, only T = 300 K) are drawn as horizontal 
arrows with respective color coding. 


as found in our simulations. In the limit R aV g —• > oo, this 
term vanishes and the MSD on planar surfaces, oc 4Z?||t, 
is restored, consistent with the vanishing non-linearity as 
observed in Fig. 2. 

By fitting Eq. (3) to our MSDs, we obtain the lateral 
diffusion constant D ||, see the inset in Fig. 2(b) for the 
fits. The results are plotted versus the Gibbs radius of 
the solutes in Fig. 3 at T = 300 K normalized by bulk dif¬ 
fusion constant Dbuik of water. For small solutes, we find 
that the parallel diffusion D \| slows down with increasing 
solute sizes which leads to parallel mobilities smaller than 
in bulk. This trend changes at a minimum at a length 
scale of R* g ss 0.45 nm beyond which the lateral diffusiv¬ 
ity monotonically rises with decreasing surface curvature 
to become bulk-like at about R* G ss 0.75 nm and finally 
saturates to a water mobility higher than in bulk, as well 
known for smooth and planar hydrophobic surfaces [35]. 
The significant length scales observed here match well 
the structural crossover length scale of ~ 0.5 nm found 
for SPC/E water estimated from solvation free energies 
of spherical model solutes [7]. 

To calculate the water diffusivity perpendicular to the 
hydrophobic surfaces, we employ the mean first-passage 
time (MFPT) analysis introduced by Hinczewski et al. 
[36, 37]. The MFPT Tf p (r,r t ) describes the mean time 
required by a molecule to travel from distance r to the so¬ 
lute to a target distance rt > r and has an exact solution 
in a Smoluchowski description in terms of the free en¬ 
ergy F(r) and diffusivity profile D±(r) [38]. Its inversion 
leads to an expression for the perpendicular diffusivity 
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Fig. 4: (a) Free energy profile F(r) for water around differ¬ 
ently sized hydrophobic solutes obtained by the Boltzmann 
inversion of the water RDFs in Fig 1. (b) MFPT curves 

Tf p (r,rt) for water molecules to reach the target distance at 
rt = 14 A given they started at a distance r < rt . (c) Perpen¬ 
dicular diffusivity profiles D±(r) obtained from Eq. (4). 


profile [36, 37] 

D±(r) = - 


3 0F(r) 


9T fp (r , r t )/dr 


f dr' e - /3F(r) (4) 

J Tmin. 


with /I -1 = ksT and r m i n being a reflective boundary 
close to the solute where F(r min ) = 10 From our 
MD simulations, we extract MFPTs T(r, r* = ro + 14 A) 
as shown in the middle panel (b) of Fig. 4 along with 
free energy profiles F(r) in the upper panel (a). The 
latter result from a simple Boltzmann inversion of the 
solute-water RDFs plus the entropic contribution F(r) = 
—ksT In g so (r) — 2fceTlnr [37]. The bottom panel (c) 
plots the resulting diffusivity profiles D±(r)/Db u ik nor¬ 
malized by bulk diffusivity for all solute sizes. The gen¬ 
eral shape of the profiles describes a maximum inside the 
first solvation shell followed by a minimum in diffusivity 
which is reached at the outside margin of the solvation 
layer, whereafter diffusivity converges towards bulk dif¬ 
fusion in an oscillatory fashion. (The final rise of the 
profiles for radii close to the target distance is an arti¬ 
fact of non-Markovian contributions [36, 37], see also the 
SI for technical details.) In our calculated profiles, an 
interesting non-monotonic curvature dependence of the 
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perpendicular water mobility near hydrophobic surfaces 
is visible in the extrema of the profiles. With growing so¬ 
lute size, the maximum of diffusivity in the first hydration 
shell continues to increase which leads to high perpen¬ 
dicular diffusivity in vicinity of weakly curved surfaces, 
whereas the minimum becomes less pronounced and al¬ 
most vanishes in the profile at the planar limit (see SI). 

We now define the perpendicular diffusion coefficient 
D±(Rmax ) in the first hydration shell by the value of the 
diffusivity curve at R max , the position of the first peak of 
the RDF. The scaled value D±(R max ) / Dbuik is presented 
next to the scaled values of D\\ in Fig. 3 versus solute size 
Rg for T = 300 K. Like the perpendicular diffusiv¬ 
ity also shows the crossover from below-to-above bulk 
mobility, but at a smaller distance of about 0.4 nm, indi¬ 
cating a preceding non-monotonic solute size dependence 
for small solutes as it is interpolated towards the limit¬ 
ing (bulk) case of a solute with vanishing radius. Since 
the magnitude of the perpendicular diffusivity slightly 
depends on the definition of what radius defines the first 
hydration shell, or if over all water molecules in this shell 
should be averaged, we compare those different defini¬ 
tions in the SI. We find similar values, in particular, all 
tested definitions rigorously reproduce the crossover from 
below-to-above bulk mobility behavior of water versus so¬ 
lute size close to the structural crossover length scale. 

Additionally Fig. 3 plots diffusivity changes upon three 
different temperatures T = 260 K, 280 K and 320 K at 
cavities at which the dynamic anomaly occurs. De¬ 
creasing temperature shifts the curves and the dynamic 
anomaly towards larger radii of curvature. Hence, the 
observed temperature dependence of diffusive dynamics 
obeys the same trend as the thermodynamic and struc¬ 
tural crossover length scale [31]. Taking into account 
entropy scaling laws for diffusion [39, 40], which scale ex¬ 
ponentially with excess entropy, suggests entropy to be 
a constitutive measure for the dynamic anomaly, thus 
firmly corroborating our conclusions. 

Ultimately, water hydrogen bond (HB) kinetics and 
mobility are intimately linked because translational dif¬ 
fusion is accompanied by processes breaking, forming 
and re-forming HBs. The cooperation of elemental dy¬ 
namic processes in water is successfully described with 
reaction-diffusion models [41-43] and was probed by 
simulations [41-44]. Here, we probe HB lifetimes r 
which are estimated by the correlation function c(t) = 
(h(t)h(0)) / (h 2 ) of the HB operator h(t), which is 1 for 
a specific pair of water molecules while they are bonded, 
and 0 otherwise [41-43]. The negative derivative k(t) = 
—dc(t)/dt is the reactive flux hydrogen bond correlation 
function. We estimate ’intermittent’ HB lifetimes by the 
zero frequency part [42], the integral of c(t), 

pOO pOO 

t = — / t ■ k(t) dt= c(t) dt . (5) 

Jo Jo 

It is significantly influenced by diffusion leading to the 



Fig. 5: Mean hydrogen bond (HB) lifetime of first-solvation 
shell water t(Rg) versus solute size scaled by the bulk value 
Tbuik (orange symbols). Average number of HBs per water 
molecule ti(Rg) scaled by its bulk value ribuik (blue symbols). 
The values corresponding to the respective measures at the 
limiting case of zero curvature (Rg —> oo, only T = 300 K) 
are drawn as horizontal arrows with respective color coding. 
The errorbars for n/nbuik were estimated from block averages 
as described in the SI. 

separation of initially bonded water molecule pairs after 
HB breaking. Without this separation, breaking a HB 
between two water molecules is reversible and the kinet¬ 
ics associated with this process are much faster than the 
irreversible rearrangements of the water hydrogen bond 
network analyzed here [42]. The correlation functions 
and details on the H-bond definition and numerical eval¬ 
uation can be found in the SI. 

The results for t(Rq) scaled by their bulk value 
Tbuik = 3.4 ps are plotted in Fig. 5. They exhibit a non¬ 
monotonic size dependence with a maximum retardation 
of about 25% versus bulk at roughly Rg — 0.45 nm, 
close to the structural crossover length scale. After 
the maximum, the retardation decreases for solute sizes 
Rg > 0.5 nm down to a remaining ss 5% near the planar 
surface. Concurrently, the average number of HBs per 
water molecule n/ribuik , also shown in Fig. 5, decreases 
monotonically. Hence, the non-monotonic translational 
mobility behavior of water is solely reflected in the ki¬ 
netic behavior of the HBs. The observed behavior recon¬ 
ciles previous, apparently contradicting results from sim¬ 
ulations [28] and theory [29, 30] on water reorientation 
times, as we find both predicted trends (increasing and 
decreasing life times), but in different solute size regimes, 
roughly separated by the important structural crossover 
scale. 

Fig. 5 also plots the T-dependence on the HB correla¬ 
tion times where a significant slowing down with decreas¬ 
ing temperature can be observed with not much change 
of the position of the maximum on the solute-size axis. 
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We note that the latter behavior cannot be related one- 
to-one to our observed diffusivity behavior in (Fig 3) due 
to the multiple timescales entering r ([42] and SI). 

In summary, we have established a firm link between 
structure and dynamics of hydration water around hy¬ 
drophobic solutes with a novel dynamic anomaly happen¬ 
ing at the well-know crossover length scale. Due to the 
fundamental importance of surface water to biomolecu- 
lar processes and function [16, 26], in particular at topo¬ 
logically heterogeneous protein surfaces [19-21, 27], our 
findings imply that nature has the means to employ local 
surface topology to mediate biological function. Hence, 
our results will help in the interpretation of experimen¬ 
tally found dynamic heterogeneities on biomolecular sur¬ 
faces [16-21]. Locally slowed down water, for instance, 
could thus fine-tune the folding kinetics of hydrophobic 
polymers and peptides [16] or may mediate the appropri¬ 
ate time scales for the association of ligands to catalyti- 
cally active sites or binding pockets [16, 22, 45]. 
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port by the Cluster of Excellence RESOLV (EXC 1069) 
funded by the Deutsche Forschungsgemeinschaft. 
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Molecular Dynamics Simulations 

The systems comprise SPC/E water as solvent and a single hydrophobic sphere fixed in the middle of the simulation 
box. The interaction between SPC/E and the model solute is mediated by a shifted Lennard Jones potential Ulj{t') = 
4e[(<T/r') 12 — (cr/r') 6 ] whereas r' = i — ro increases the solute size. The parameters e = 1.197 kJ/mol and a = 3.768 A 
are those of a united-atom methane molecule as in Ref. [1]. Ten seperate simulations are performed with shift radii 
ro = OA, 1 A, 2 A, 4 A, 6 A, 8 A, 10 A, 12.5 A, 15 A and 17.5 A. For shift radii ro < 6 A the simulation box containes 
6000 SPC/E water molecules and for simulations with ro > 6 A the number of solute molecules is doubled to 12000 
to avoid finite size effects. 

To model the corresponding hydrophobic surface without curvature the walls in the x-y-plane of a simulation box 
of size 5.67 x 5.67 x 6.54 nm 3 are set to interact with water with a 12-6 potential in the z-direction Ui 2 -s(z) = 
4e[(cr/z) 12 — {o/z) 6 ] with same interaction parameters e = 1.197kJ/mol and a = 3.768A. 

Each system is equilibrated in NPT ensemble for 100 ps at ambient conditions, namely T = 300 K and P = 1 bar. 
After equilibration the box lengths of the cubic boxes are roughly 5.65 nm and 7.15 nm for the small (6000 SPC/E) 
and big (12000 SPC/E) setups respectively. In order to gather enough statistics to probe long-time dynamics within 
the thin solvation shell comprised of 0(10) to 0(1O 2 ) water molecules subsequent production runs have a length 
of 200 ns for the small systems and 100 ns for the big systems with a time step of 2fs and writing period of 200 fs. 
Energy divergence is prevented by simulation in NVT ensemble (T = 300 K) applying the Nose-Hoover thermostat 
with a coupling period of 1 ps. The Nose-Hoover is implemented by an extra degree of freedom introduced as a heat 
bath within the hamiltionian of the simulation and thus has the advantage of creating physically more realistic NVT 
ensembles [5, 6]. 

Additional simulations for the cavities with ro = 0A, lA, 2 A and 4 A were performed in order to resolve the 
temperature dependence of our observations. Production runs with T = 260 K, 280 K and 320 K were performed up 
to 100 ns. 

All simulations are performed with the GROMACS 4.5.4 package [4]. 


Planar Limit 

The spatial distribution function in figure 1 of water at planar walls with the interaction parameters mentioned 
above does not exhbit a clear first and second minimum. The missing extrema hinder a clear definition of a first and 
second solvation shell as it can be observed from the RDFs (main text Fig. 1) for curved surfaces. 

One further simulation of a wall interacting with a doubled interaction strength e = 2.394 kJ/mol was conducted 
to generate clear extrema (Fig. 1) and hence definitions for the solvation layers. Clearly, the found borders to the 
solvation shells only approximate values that correspond to values that would consistently define hydration layers for 
the original interaction parameters. 

Another notable difference to simulations of water at curved surfaces arrises from the lack of long-range correction 
for the interaction in the planar geometry. This creates smaller water densities and larger mobilities than in bulk [3]. 

Both discrepancies introduce minor consistency deviations between the simulations with the curved surfaces and 
the planar limit but never conflict with our interpretation. 
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Fig. 1: Normalized water distribution at hydrophobic walls. 




Fig. 2: Left panel: Spatially resolved perpendicular diffusion profile of water at the hydrophobic wall. Right panel: (red, 
solid) Average perpendicualr diffusion constant (D±(r)) (eq. 1) inside the solvation layer of differently sized hydrophobic model 
solutes. The red horizontally dashed line is the corresponding value at the planar limit (Rg —> oo) for (D±(r)) (green, dashed) 
Excess perpendicular diffusion constant D exc + 1 (eq. 2) quantifying an average disruption of the diffusion of water near 
hydrophobic solutes. The reference value of bulk diffusivity is drawn as grey horizontally dashed line as isoline with value 1. 


MFPT method and perpendicular diffusivity 

In the main text we briefly describe the MFPT method from Hinczewski et al. [2] which takes the free energy profile 
F[r ) and the spatial derivative of the MFPT curve Tf p (r,r t ). Both are obtained equivalently to the application of 
pair diffusion of water and methane by Sedlmeier et al. [3]. The integral J’ dre - ^^ is evaluated numerically 
starting from the reflective boundary r m i n , which we define by the largest distance from the respective solute where 
F(r m i n ) > 10 ks T. The derivative dTf p (r,rt)/dr is determined by fitting a linear function to the function values 
around r — Sr < r < r + Sr whereas 5r = 0.05 nm. The closest distance r evaluated for the derivative, and hence 
diffusivity profile value closest to the solute, is r S Rg where statistics of Tf p (r,r t ) allow a reasonable fit in the 
mentioned range of Sr. 

The left panel of Fig.2 plots the diffusivity profile of water motion perpendicular to the hydrophobic planar model 
surface. The diffusivity right at the surface decreases rapidly from high mobility values to a global minimum proceeding 
into a maximum which is located within the hydration water layer. With growing distance to the wall the diffusivity 
then decreases to bulk diffusion. The minimum at the outer margin of the solvation shell which is visible in each 
profile in Fig.4 of the main text almost fully vanishes for the profile of the planar surface here. 

The right panel in Fig.2 shows two further measures quantifying the perpendicular diffusivity in vicinity of the 
hydrophobic model solutes. The red curve is an average diffusivity inside the first solvation shell 

(D±{r)) = f drg so (r)-D ± (r). (1) 

JRma* 
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Fig. 3: Perpendicular diffusion at cavities with ro = 0 A, 1 A, 2 A and 4 A with temperatures T = 260 K (upper left), T = 280 K 
(upper right), T = 300 K (bottom left) and T = 320 K (bottom right). 


It captures the non-monotonicity from the diffusivity profiles slightly decreasing below bulk diffusion near small solutes 
and subsequently increases converging towards the value for perpendicular diffusion at the planar interface. 

We also calculate an ’’excess” perpendicular diffusivity D exc which quantifies an average change of the diffusivity 
constant in the solute’s surrounding water 

/»00 

D eX c= / drr 2 g so (r)-[D ± (r)/D bulk -l] . (2) 

■JRmau: 

The plot shows 1 +D exc {Rc)- Certainly the non-monotonic curvature dependence in the diffusivity profiles is smoothly 
captured in this quantity aswell. In general an excess diffusivity constant can be of particular interest in order to 
estimate water diffusivity changes in dilute solutions of hydrophobic solutes. 


Hydrogen Bond Time Correlation 

We analyse the autocorrelation function c(t ) of the H-bond existence inside the first hydration layer using an in- 
house analysis code. It uses a geometric H-bond definition with a minimum donor-acceptor distance duB < 3.5 A and 
donor-hydrogen-acceptor angle 9hb > 150°. 

Fig. 4 plots the correlation functions c(f) for hydrogen bonds inside the first hydration layer. Hydrogen bonds were 
selected at time 0 if they were intact and both, acceptor and donor molecules, resided in the first hydration shell. 
We evaluate c(t) up to correlation times of 100 ps until when it has decayed to magnitudes of 10 -3 . Therefore, the 
integral over the correlation function in this time window provides a reasonable lower boundary to the actual HB 
lifetimes. 

In addition the average number of hydrogen bonds per water molecule n within the solvation layer was counted and 
presented in the main text. The error of the average number of hydrogen bonds per molecule is estimated by block 
averaging. Using a frequency of a 1000 ps the average number of hydrogen bonds within a block of lOOps length was 
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Fig. 4: Correlation function c(t) of the hydrogen bonding operator h(t) of hydrogen bonds on water molecules inside the first 
solvation shell with temperatures T = 260K (upper left), T = 280K (upper right), T = 300K (bottom left) and T = 320K 
(bottom right). Only at T — 300K all cavity sizes were evalutated. For the rest the analysis was limited to cavities with 
ro = 0 A, 1 A, 2 A and 4 A. 


calculated, giving 100 uncorrelated block averages nBlock- These were used to approximate the standard deviation 
error by Sn = [T, Biocks ( n - n B i oc k) 2 / 100] 1 / 2 . 


Temperature dependence of hydrogen bond dynamics 

The hydrogen bond correlation function c(t) contains all three time scales of the reaction-diffusion model from 
Luzar and Chandler [7] given by 


d 

—p(r, t) = D\/ 2 p(f ., t) + 5(f) kc(t ) — 5(f) k'n(t ) 
ot 

= Z?V 2 / o(r, t ) + k(t)5r, 


where the rates k and k! are those for breaking and reformation of a hydrogen bond, respectively. n(t) is proportional 
to the diffusivity since n(t) oc p(0,t) [8] (and p(r,t) obeying Fick’s law). 

Thus, as we used the reactive flux k(t) = —c(t) to evaluate the zero frequency part 


T 



' + k'n(t ) 
k 


dt, 


the correlation time r contains the factor Dfc'/fc, which precludes a direct one-to-one connection to the temperature 
dependence of the diffusivity. Certainly a direct relation of temperature dependence between hydrogen bond dynamics 
and diffusivity is restored if the time scales were extracted separately, which, however, is not a straightforwared task. 
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